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Abstract 

Quantum computing shows great promise for the solution of many dif- 
ficult problems, such as the simulation of quantum systems and the fac- 
torisation of large numbers. While the theory of quantum computing is 
well understood, it has proved difficult to implement quantum computers 
in real physical systems. It has recently been shown that nuclear magnetic 
resonance (NMR) can be used to implement small quantum computers us- 
ing the spin states of nuclei in carefully chosen small molecules. Here we 
demonstrate the use of an NMR quantum computer based on the pyrim- 
idine base cytosine, and the implementation of a quantum algorithm to 
solve Deutsch's problem (distinguishing between constant and balanced 
functions). This is the first successful implementation of a quantum algo- 
rithm on any physical system. 
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I. INTRODUCTION 



In 1982 Feynman pointed out that it appears to be impossible to efficiently sim- 
ulate the behaviour of a quantum mechanical system with a computed. This 
problem arises because the quantum system is not confined to its eigenstates, 
but can exist in any superposition of them, and so the space needed to describe 
the system is very large. To take a simple example, a system comprising N 
two level sub-systems, such as N spin-| particles, inhabits a Hilbert space of 
dimension 2^, and evolves under a series of transformations described by ma- 
trices containing 4^ elements. For this reason it is impractical to simulate the 
behaviour of spin systems containing more than about a dozen spins. 

The difficulty of simulating quantum systems using classical computers sug- 
gests that quantum systems have an information processing capability much 
greater than that of corresponding classical systems. Thus it might be possi- 
ble to build quantum mechanical computera^EIEEl which utilise this information 
processing capability in an effective way to achieve a computing power well be- 
yond that of a classical computer. Such a quantum conmuter could be used to 
efficiently simulate other quantum mechanical systema^^, or to solve conven- 
tional mathematical problemJ^ which suffer from a similar exponential growth in 
complexity, such as factoring^. 

Considerable progress in this direction has been made in recent years. The ba- 
sic logic elements necessary to carry out quantum computing are well understood, 

and quantum algorithms have been developed, both for simple demonstration 
problemfP^E 

and for more substantial problems such as factoring^^. Exper- 
imental implementation of a quantum computer has, however, proved difficult. 
Much effort has been directed towards implementing quantum computers using 
ions trapped by electric and magnetic fields- , and while this approach has shown 
some success^, it has proved difficult to progress beyond computers containing 
a single two-level system (corresponding to a single quantum bit, or qubit). 

Recently two separate approaches have been describecP*^^ for the imple- 
mentation of a quantum computer using nuclear magnetic resonanc tP^l (NMR). 
These approaches shows great promise as it has proved relatively simple to inves- 
tigate quantum systems containing two or three Here we describe 
our implementation of a simple quantum algorithm for solving Deutsch's prob- 
lemj^^El on a two qubit NMR quantum computer. 



II. QUANTUM COMPUTERS 

All current implementations of quantum computers are built up from a small 
number of basic elements. The first of these is the qubit, which plays the same 
role as that of the bit in a classical computer. A classical bit can be in one of 
two states, or 1, and similarly a qubit can be represented by any two level 
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system with eigenstates labelled |0) and |1). One obvious implementation is to 
use the two Zeeman levels of a spin-| particle in a magnetic field, and we shall 
assume this implementation throughout the rest of this paper. Unlike a bit, 
however, a qubit is not confined to these two eigenstates, but can in general exist 
in some superposition of the two states. It is this ability to exist in superpositions 
which makes quantum systems so difficult to simulate and which gives quantum 
computers their power. 

The second requirement is a set of logic gates, corresponding to gates such as 
and, OR and not in conventional computers!^ Quantum gates differ from their 
classical counterparts in one very important way: they must be reversible^"^*^. 
This is because the evolution of any quantum system can be described by a 
series of unitary transformations, which are themselves reversible. This need for 
reversibility has many consequences for the design of quantum gates. Clearly for 
a gate to be reversible it must be possible to reconstruct the input bits knowing 
only the design of the gate and the output bits, and so every input bit must 
be in some sense preserved in the outputs. One trivial consequence of this is 
that the gate must have exactly as many outputs as inputs. For this reason it is 
obvious that gates such as and and OR are not reversible. It is, however, possible 
to construct reversible equivalents of and and OR, in which the input bits are 
preserved. 

Just as it can be shown that one or more gates (such as the nand gate) are 
universal for classical computing^! (that is, any classical gate can be constructed 
using only wires and nand gates), it can be shown that certain gates or com- 
binations of gates are universal for quantum computing. In particular it can be 
shown^l that the combination of a general single qubit rotation with the two bit 
"controlled-NOT" gate (cnot) is universal. Furthermore it is possible to build a 
reversible equivalent of the NAND gate, and thus to implement any classical logic 
operation using reversible logic. 

Single qubit rotations are easily implemented in NMR as they correspond to 
rotations within the subspace corresponding to a single spin, and such rotations 
can be achieved using radio frequency (RF) fields. One particularly important 
single bit gate is the Hadamard gate, which performs the rotational transforma- 
tion 

10) 
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The Hadamard operator can thus be used to convert eigenstates into superpo- 
sitions of states. Similarly, as the Hadamard is self-inverse, it can be used to 
convert superpositions of states back into eigenstates for later analysis. 

Two-bit gates correspond to rotations within subspaces corresponding to two 
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input output 
| I 
10 1 

10 11 

11 10 



Table 1: The truth table for the cnot gate. The first qubit (the control qubit) 
is unchanged by the gate, while the second qubit is flipped if the control qubit is 
in state 1, effectively implementing an XOR gate. 

spins, and thus require some kind of spin-spin interaction for their implementa- 
tion. In NMR the scalar spin-spin coupling (J-coupling) has the correct form, 
and is ideally suited for the construction of controlled gates, such as cnot. This 
gate operates to invert the value of one qubit when another qubit (the control 
qubit) has some specified value, usually |1); its truth table is shown in tabled 

Finally it is necessary to have some way of reading out information about the 
final quantum state of the system, and thus obtaining the result of the calcula- 
tion. In most implementations of quantum computers this process is equivalent 
to determining which of two eigenstates a two level system is in, but this is not 
a practical approach in NMR. It is, however, possible to obtain equivalent infor- 
mation by exciting the spin system and observing the resulting NMR spectrum. 
Different qubits correspond to different spins, and thus give rise to signals at dif- 
ferent resonance frequencies, while the eigenstate of a spin before the excitation 
can be determined from the relative phase (absorption or emission) of the NMR 
signals. 

III. THE DEUTSCH ALGORITHM 

Deutsch's problem in its simplest form concerns the analysis of single bit binary 
functions: 

f{x) :B^B, (2) 

where B = {0, 1} is the set of possible values for a single bit. Such functions 
take a single bit as input, and return a single bit as their result. Clearly there 
are exactly four such functions, which may be described by their truth tables, 
as shown in table El These four functions can be divided into two groups: the 
two "constant" functions, for which f(x) is independent of x (f 00 and /n), and 
the two "balanced" functions, for which f(x) is zero for one value of x and unity 
for the other (/ i arid /io). Given some unknown function / (known to be one 
of these four functions), it is possible to determine which of the four functions 
it is by applying / to two known inputs: and 1. This procedure also provides 



4 



X 


/oo(aO 






h\{x) 











1 


1 


1 





1 





1 



Table 2: The four possible binary functions mapping one bit to another. 

enough information to determine whether the function is constant or balanced. 
However knowing whether the function is constant or balanced corresponds to 
only one bit of information, and so it might be possible to answer this question 
using only one evaluation of the function /. Equivalently, it might be possible to 
determine the value of /(0) © /(l) using only one evaluation of /. (The symbol 
© indicates addition modulo 2, and for two one bit numbers, a and b, a@b equals 
if a and b are the same, and 1 if they are different.) In fact this can be achieved 
as long as the calculation is performed using a quantum computer rather than a 
classical one. 

Quantum computers of necessity use reversible logic, and so it is not possible 
to implement the binary function / directly. It is, however, possible to design a 
propagator, Uf, which captures / within a reversible transformation by using a 
system with two input qubits and two output qubits as follows 

\x)\y)-%\x)\y®f(x)). (3) 

The two input bits are preserved (x is preserved directly, while y is preserved by 
combining it with f(x), the desired result), and so Uf corresponds to a reversible 
transformation. Note that for any one bit number a, © a = a, and so values 
of f(x) can be determined by setting the second input bit to 0. Using this 
propagator and appropriate input states it is possible to evaluate /(0) and /(l) 
using 

|0)|0>-^|0>|/(0)> (4) 

and 

|1>|0)^|1)|/(1)). (5) 

The approach outlined above, in which the state of a quantum computer is de- 
scribed explicitly, swiftly becomes unwieldy, and it is useful to use more compact 
notations. One particularly simple approach is to use quantum circuits,^! which 
may be drawn by analogy with classical electronic circuits. In this approach lines 
are used to represent "wires" down which qubits "flow", while boxes represent 
quantum gates which perform appropriate unitary transformations. For example, 
the analysis of / can be summarised by the circuit shown in figure ^ 
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Figure 1: Quantum circuit for the analysis of a binary function /. 



So far this is simply using a quantum computer to simulate a classical com- 
puter implementing classical algorithms. With a quantum computer, however, it 
is not necessary to start with the system in some eigenstate; instead it is possible 
to begin with a superposition of states. Suppose the calculation begins with the 
second qubit in the superposition (|0) — |l))/v2. Then 
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|0 ) - ]i) 
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(We have used the fact that © a = a, as before, while l©a = lifa = 
and if a = 1). The value of f(x) is now encoded in the overall phase of the 
result, with the qubits left otherwise unchanged. While this is not particularly 
useful, suppose the calculation begins with the first qubit also in a superposition 
of states, namely (|0) + \ l))/y/2. TherP 
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with the first qubit ending up in the superposition (|0)±| 1))/ v2, with the desired 
answer (/(0) ©/(l)) encoded as the relative phase of the two states contributing 
to the superposition. This relative phase can be measured, and so the value of 
/(0) © /(l) (that is, whether / is constant or balanced) has been determined 
using only one application of the propagator Uf, that is only one evaluation of 
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Figure 2: A quantum circuit for solving Deutsch's problem. 



the function /. 

This approach can be easily implemented using a quantum circuit, as shown 
in figure |21 This circuit starts off from appropriate eigenstates, uses Hadamard 
transformations to convert these into superpositions, applies the propagator Uf 
to these superpositions, and finally use another pair of Hadamard transforms to 
convert the superpositions back into eigenstates which encode the desired result. 



IV. IMPLEMENTING THE DEUTSCH ALGORITHM IN NMR 



The Deutsch algorithm can be implemented on a quantum computer with two 
qubits, such as an NMR quantum computer based on two coupled spins. First 
it is necessary to show that the individual components of the quantum circuit 
can be built. It is convenient to begin by writing down the necessary states and 
operators using the product operator basis seiP^^ normally used in describing 
NMR experiments (this basis is formed by taking outer products between Pauli 
matrices describing the individual spins, together with the scaled unit matrix, 
\E). 

The initial state, |"0oi) = can be written as a vector in Hilbert space, 



|-0O1> 
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(8) 



but this description is not really appropriate. Unlike other implementations, 
an NMR quantum computer comprises not just a single set of spins but rather 
an ensemble of spins in a statistical mixture of states. Such a system is most 
conveniently treated using a density matrix, which can describe either a mixture 
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or a pure state; for example 



P01 = |^O1><-0O1 



/0 0\ 

10 



\0 0/ 



(9) 



This density matrix can be decomposed in the product operator basis as p i — 
(I z — S z — 2/ 2 S' 2 + \E)/2. Ignoring multiples of the unit matrix (which give rise 
to no observable effects in any NMR experiment) this can be reached from the 
thermal equilibrium density matrix (I z + S z ) by a series of RF and field gradient 
pulsed^. 

The unitary transformation matrix corresponding to the Hadamard operator 
on a single spin can be written as 



This corresponds to a 180° rotation around an axis tilted at 45° between the z and 
i ... ff^l 

Such a rotation can be achieved directly using an off resonance pulse, 1 **"" 1 

or using a three pulse sandwicfP*^* such as 45° — 180°— 45°_ y . Even more simply the 

Hadamard can be approximated by a 90° pulse. While this is clearly not a true 

Hadamard operator (for example, it is not self inverse), its behaviour is similar 

and it can be used in some cases: for example, it is possible to replace the first 

pair of Hadamard gates in the circuit for the Deutsch Algorithm (figure EJ by 90° 

pulses and the second pair of gates by 90°_ y pulses. Clearly it is possible to apply 

the Hadamard operator either to just one of the two spins (using selective soft 

RF pulses^) or to both spins simultaneously (using non-selective hard pulses). 

The unitary transformations corresponding to the four possible propagators 

Uf are also easily derived. Each propagator corresponds to flipping the state of 

the second qubit under certain conditions as follows: U 00 , never flip the second 

qubit; Uqi, flip the second qubit when the first qubit is in state one; Uio, flip 

the second qubit when the first qubit is in state zero; Un, always flip the second 

qubit. The first and last cases are particularly simple, as Uqo corresponds to doing 

nothing (the identity operation), while Un corresponds to inverting the second 

spin (a conventional not gate, or, equivalently, a 180° pulse). The second and 

third propagators correspond to controlled-NOT gates, which can be implemented 

using spin-spin couplings. For example Uqi is described by the matrix 




(10) 
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(11) 
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which can be achieved using the pulse sequence 

90 S y - couple - 90 I z - 90 S_ 2 - 90 S_ y (12) 

where 90 S y indicates a 90° pulse on the second spin, couple indicates evolution 
under the scalar coupling Hamiltonian, ttJjs2I z S z , for a time 1/2 Jjs, and 90 I z 
and 90 S- z indicate either periods of free precession under Zeeman Hamiltonians 
or the application of composite z-pulsesP^EI] Similarly Uiq can be achieved using 
the pulse sequence 

90 S y - couple - 90 I z - 90 S z - 90 (13) 

The pulse sequences described above can be implemented in many different 
ways, as different composite ^-pulses can be used, the order of some of the pulses 
can be varied, and in some cases different pulses can be combined together. We 
chose to use the implementation 

90 S y - 1/4 J/s - 180 x . - 1/4 J/s - 180 x - 90 I v - 90 I x - 90_„ - 90 S± x (14) 

where pulses not marked as either I or S were applied to both nuclei. The phase 
of the final pulse distinguishes Uoi (for which the final pulse is S +x ) from Uio (for 
which it is S- x ). 

Finally it is necessary to consider analysis of the final state, which could in 
general be one of the four states p o, p i> Pith or Pn- m order to distinguish 
these states it is necessary to apply a 90° pulse and observe the NMR spectrum. 
The final NMR signal observed from spin I is I x if the spin is in state 0, and 
— I x if it is in state 1. For a computer implementing the Deutsch algorithm the 
final detection 90° pulses cancel out the two final pseudo-Hadamard 90°.^ pulses, 
and thus all four pulses can be omitted (see figure The final NMR signal 
observed is either \l x — \S X (corresponding to p i) or ~~ — ^S x (corresponding 
to pn). Hence it is simple to determine the value of /(0)©/(l) (that is, determine 
whether the function is constant or balanced) by determining the relative phase 
of the signals from the two spins. 

V. EXPERIMENTAL RESULTS 

In order to demonstrate the results described above, we have constructed an NMR 
quantum computer capable of implementing the Deutsch algorithm. For our two- 
spin system we chose to use a 50 mM solution of the pyrimidine base cytosine 
in D2O; rapid exchange of the two amine protons and the single amide proton 
with the deuterated solvent leaves two remaining protons forming an isolated 
two-spin system. All NMR experiments were conducted at 20°C and pH* = 7 
on a home-built NMR spectrometer at the Oxford Centre for Molecular Sciences, 



9 



|0> 
10} 









9o; 














90; 











+x 



(b) 
10) 

|1) 





90°, 




Uf 














90°, 











—x 



Figure 3: Modified quantum circuits for the analysis of binary functions on an 
NMR quantum computer, (a) A circuit for the classical analysis of /(0); the 
normal circuit (see figured]) is followed by 90° pulses to excite the NMR spectrum. 
Clearly f(l) can be obtained in a similar way. (b) A circuit for the implementation 
of the Deutsch algorithm, with Hadamard operators replaced by 90°^ pulses. The 
final 90° excitation pulses cancel out the 90°_^ pulses, and thus all four pulses can 
be omitted. 
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Figure 4: Experimental implementation of an algorithm to determine /(0) on an 
NMR quantum computer, (a) The result of applying Uf 00 ] as this propagator is 
the identity matrix this spectrum can also serve as a reference. The lefthand pair 
of signals corresponds to the first spin (J), while the pair on the right hand side 
correspond to the second spin (S) . Note that the signals from both spins (which 
are in state |0), the ground state) are in absorption, (b) The result of applying 
Uf 01 ] both sets of signals are still in absorption, as /(0) =0 for this function, (c) 
The result of applying Uf 10 ; the signals from spin S are now in emission, since 
/(0) = 1 for this function, (d) The result of applying Uf u ; the signals from spin 
S are once again in emission as expected. 

with a 1 H operating frequency of 500 MHz. The observed J-coupling between the 
two protons was 7.2 Hz, while the difference in resonance frequencies was 763 Hz. 
Selective excitation was achieved using GaussiarP^l so ft pulses incorporating a 
phase ram to allow excitation away from the transmitter frequency. During 

a selective pulse the other (unexcited) spin continues to experience the main 
Zeeman interaction, resulting in a rotation around the z-axis, but the length of 
the selective pulses can be chosen such that the net rotation experienced by the 
other spin is zero. The residual HOD resonance was suppressed by low power 
saturation during the relaxation delays. 

This system can be used both for the implementation of classical algorithms 
to analyse /(0) and /(l) and for the implementation of the Deutsch algorithm; as 
shown in figure|3]the pulse sequences differ only in the placement of the 90° pulses. 
The results for the classical algorithm to determine /(0) are shown in figure |U 
The lefthand pair of signals corresponds to the first spin (J), while the pair on the 
right hand side correspond to the second spin (S); the (barely visible) splitting 
in each pair arises from the scalar coupling Jj$. In this experiment the value 
of /(0) is determined by setting both spins I and S into state |0), performing 
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(a) 



(b) 



(c) 
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(d) 



Figure 5: Experimental implementation of an algorithm to determine /(l) on an 
NMR quantum computer; in this case the algorithm starts with spin / in the 
excited state, |1), and so signals from spin I are in emission. For details of the 
labelling see figure |H 



the calculation, and then measuring the final state of spin S; spin / should 
not be affected, and so should remain in state |0). The phase of the reference 
spectrum (a) was adjusted so that signals from spin / appear in absorption, and 
the same phase correction was applied to the other three spectra. The state of 
a spin after a calculation can then be determined by determining whether the 
corresponding signals in the spectrum are in absorption (state |0)) or emission 
(state |1)). As expected spin / does indeed remain in state |0), while the value 
of /(0) (determined from spin S) is for Uf 00 and Uf 01 but 1 for Uf 10 and Uf u . 

Clearly our NMR quantum computer is capable of implementing this classical 
algorithm, as it is simple to determine /(0). The other value, /(l), can be 
determined in a very similar way (see figure EJ). In this case spin / remains in 
state |1), while /(l) equals for Uf 00 and Uf 10 and equals 1 for Uf m and Uf u . 
There are, however, several imperfections visible in the results. 

First, the signals are not perfectly phased: rather than exhibiting pure absorp- 
tion or pure emission lineshapes the signals have more complex shapes, including 
dispersive components. These arise from the difficulty of implementing perfect se- 
lective pulses, which effect the desired rotation at one spin while leaving the other 
spin entirely unaffected. Similarly the selective pulses will not perfectly suppress 
J-couplings during the excitation, leading to the appearance of antiphase contri- 
butions to the lineshape. Any practical selective pulse will be imperfect, and so 
will result in systematic distortions in the final result. Note that these distortions 
are most severe in cases (b) and (c), where the propagator is complex, containing 
a large number of selective pulses. Interestingly the distortions are also more 
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(a) 
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(c) 



(d) 



Figure 6: Experimental implementation of a quantum algorithm to determine 
/(0) © /(l) on an NMR quantum computer. In this case the result can be read 
out on spin /, that is using the signals on the left of the spectrum. For details of 
the labelling see figure 0] As expected the I spin is inverted when the function 
is balanced (/oi or /i ), but not when the function is constant (/oo or fu)- 



severe for the measurement of /(0) (figureHJ) than for /(l) (figure^; there is no 
simple explanation for this effect, which is due to the complex interplay of many 
selective pulses. We are currently seeking ways to minimise these effects. 

Secondly, the signal intensities vary in different cases; as before the signal 
loss is most severe in cases (b) and (c), corresponding to complex propagators. 
This is in part a consequence of imperfect selective pulses, as discussed above, 
but may also indicate the effects of spin relaxation, that is decoherence of the 
states involved in the calculation. Decoherence is a fundament al p roblem . and 
may ultimately limit the size of practical quantum computers 12511213123 although 
a variety of error correction techniqueP^^SEDl have been devised to overcome 
it. 

These imperfections are not a major problem in our NMR quantum computer, 
as it is still easy to determine the state of a spin. However our computer is small, 
and the programs run on it are short (that is, they contain a small number of 
logic gates); if more complex programs are to be run on larger computers then 
these imperfections must be addressed. 

The results of implementing the Deutsch quantum algorithm are shown in 
figure EJ In this case the result (/(0) © /(l)) can be read from the final state of 
spin /, while spin S remains in state |1). As expected spin / is in state |0) for 
the two constant functions (/oo an d fu), but in state |1) for the two balanced 
functions (foi and /io). Once again a number of imperfections are visible, though 
in this case they appear to be most severe in the case of Uf n . 
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VI. SUMMARY 



We have demonstrated that the isolated pair of 1 H nuclei in partially deuterated 
cytosine can be used to implement a two qubit NMR quantum computer. This 
computer can be used to run both classical algorithms and quantum algorithms, 
such as that for solving Deutsch's problem (distinguishing between constant and 
balanced functions). This is the first successful implementation of a quantum 
algorithm on any physical syste nP. 

This result confirms that NMR shows great promise as a technology for the 
implementation of small quantum computers. Difficulties do exist, largely as a 
result of the large number of selective pulses involved in the implementation of 
quantum gates, but we are currently seeking ways to overcome these problems. 
Even with the current level of errors it should be possible to build a three qubit 
computer capable of implementing more complex logic gates and algorithms. 
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